Approach to the extremal limit of the Schwarzschild— de Sitter black hole 
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The quasinormal-mode spectrum of the Schwarzschild-de Sitter black hole is studied in the limit 
of near-equal black-hole and cosmological radii. It is found that the mode frequencies agree with the 
Poschl-Teller approximation to one more order than previously realized, even though the effective 
potential does not. Whether the spectrum approaches the limiting one uniformly in the mode index 
is seen to depend on the chosen units (to the order investigated). A perturbation framework is set 
up, in which these issues can be studied to higher order in future. 
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Since Schwarzschild's original black-hole solution of 
the Einstein equations, several extensions have been 
found, including those with nonzero rotation or elec- 
tric charge [j , or the Schwarzschild-de Sitter (SdS) hole 
with a finite outer 'cosmological' horizon Q. Each of 
these leads to a corresponding generalization of the wave 
equation governing the linearized perturbations Q ■ Re- 
laxing the requirement of physical significance, one can 
also generalize this wave equation by staying with the 
Schwarzschild spacetime, but taking the spin s of the 
perturbing field to be formally continuous (gravity waves 
correspond to s = 2). In all these cases, one faces 
the question how the Schwarzschild quasinormal-mode 
(QNM) spectrum is modified. This question is compli- 
cated by said spectrum being infinite, with two strings 
of high-damping modes — currently undergoing intense 
scrutiny £| — extending to uj — —ioo; therefore, one can 
have nontrivial interplay between the Schwarzschild and 
large-mode-index limits p|. 

Clearly, the study of these and analogous situations 
yearns for a solvable zeroth-order model. Recently, one 
has become available by the demonstration Q that the 
SdS wave equations reduce to the familiar Poschl-Teller 
(PT) model 7] in the opposite 'extremal' limit. While 
the PT potential was already being used as a toy model 
for gravitational problems y|, the derivation puts this 
on a firm footing and in particular provides a class of 
physically motivated perturbations, for which moreover 
numerical data are already available In this report, 
these perturbing potentials are derived in detail, and a 
beginning is made with their study. 

The SdS spacetime is characterized by a black-hole 
horizon, set to unity here, and a cosmological horizon c, 
the only nontrivial parameter. For each angular momen- 
tum £, the perturbations obey a wave equation 
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where the tortoise coordinate is given by dr* — dr/f(r) 
in terms of the warp factor 



/ = 



(r - l)(c - r)(r + 1 + c) 
rA 
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tion (i.e., on s), the effective potential reads 
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with M = c(l+c)/2A In Ref. 9], the ensuing QNM 
spectra are compared to a PT approximation for all c. 
However, the used fit amounts to reproducing the fun- 
damental (complex) frequency to high accuracy, while 
in fact it is the high-order modes that — at least in the 
Schwarzschild limit — qualitatively have PT character. 
Thus, there is room for further analytical study of the 
PT approximation. 

For c — ► oo, O - © readily reduce to the familiar 
Regge-Wheeler equation for the Schwarzschild space- 
time. Here, however, we are interested in the limit c [ 1, 
where the inner and outer radii merge. Let us scale 
c = 1 + e, r = 1 + e£ (0 < £ < 1), r* = x/e, V = e 2 W, 
and u> = ev. The new tortoise coordinate x = e 2 fd£/ f is 



(3+3e+e 2 ) 
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(3+e)(3+2e) 



ln(3+e+e£) 



(4) 



In the scaled variables, expansion in e is straightforward. 
Inverting x(£), one finds 
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with £o = l/(l+e~ x ). Although we intend to carry our 
expansions to second order, we will only do so after one 
further simplification. Namely, consider electromagnetic 
perturbations. Inserting JSJ) into |JS}, one obtains 



l-2e 



£(£+!) 4ch 2 (:r/2) 
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with A = 1 + c + c 2 . Depending on the type of perturba- 



sh(x/2) 

ch 3 (x/2) ' 
(6) 

One observes that most of the (D(e)-corrections can in fact 
be absorbed into a further redefinition of the coordinate. 
In other words, it is — unsurprisingly — advantageous to 
fit V c \ to a PT potential with e-dependent position and 
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width. Thus, we set x = (l+e)y 



le,W 



U/(l+e) 2 , and v = /i/(l+e), leading to [6 = l/(l+e~ y )} 



^ 2e ln(2ch|) + 2 £ 2 2(l-e^)ln 2 (2ch|) + (9-3e y ) In (2 ch|) - | (l+e y )y + 16 + 6e y 
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+ 0(e 3 ), (7) 



4ch 2 (y/2) 6 ch 3 (y/2) 



Z(ef-4+e-y)ln 2 (2ch|) + (4-5e y +e- y ) ln^2 ch|) + 5sh(y)y 
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Since the 0(e)-term is now odd m these coordinates there are no O(e)-corrections at all. However, this C(e)-term 
in U e \ can contribute in 0(e 2 ), see the Appendix. For studying perturbation to the latter order, the 0(e 2 )-term in 
U e \ can be safely symmetrized, so that the numerator becomes 



8ch 2 (!)-12] ln 2 (2ch|) - 8sh 2 (f ) ln( 2 ch|) + 10ch(|) sh(|)y - 20ch 2 (f ) + 15 . 
For scalar and gravitational perturbations, the analogous potentials read 



u„ = U R 



e s%/2) 
4ch 3 (2//2) 



4ch 2 (y/2) 
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£ 2 A(A-e v -e- y ) \n\2 ch| J + (16e H -32+4e- J/ ) In (2 ch-J - 10 sh(y)y + 17e^ - 13 + 6e~ v 



144 



ch 4 (y/2) 



+ 0(e 3 ). (11) 



Again, the leading corrections are odd. The esh(y/2)x 
ch~ (y/2)-terms could be eliminated also here by making 
the x t— ► y shift s- and ^-dependent, but this will not 
be pursued. Only after the first-order terms in IjlOfl - 
(|llfl will thus have been maximally simplified does it look 
meaningful to symmetrize the second-order terms. 
Defining 



9 = 




(12) 



in all cases one has [i n = |[g — (n+|)i] + 0(e 2 ), n = 
0,1,2,..., so that w„ = i(e-e 2 )[(7-(n+i)i] + 0(e 3 ). An 
important quantity in SdS is the surface gravity |(| 
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Since < e < 00 but < K < \, it is heuristically 
advantageous to reparametrize the problem in terms of k, 
inverting (O as e = 2k + |k 2 + 4r« 3 H , so that 



w n = «(l-|«)[g - (n+|)t] + 0(k 3 ) 



(14) 



The factor 1 — |/c is also the leading one in (19) of 
but here the general ^-dependence has been obtained ex- 
actly. For a merely moderately small k = 0.01411708, 
JHJ yields 



/ 0.0270803 - 0.006992Hz, 
w °' gr w \0.0436657 - 0.0069921H , 

while numerically one has [9j 



£ = 2 
^ = 3 



0.0270837 - 0.00699266i , £ = 2 
0.0436710 - 0.00699280i , £ = 3 



(15) 



(16) 



For both values of £, the error has been verified to de- 
crease cubically with n. 

Thus, uj n /e does have O(en)-corrections [u> n /K has 
C(Kn)-corrections]. In these units one could say that 
there is 'non- uniform convergence to the limit', while in 
fact /j. n does not have any corrections at all in this or- 
der. Of course, it is possible (even likely) that the 0(e 2 )- 
corrections to ji n are unbounded in n. The above, plus 
the general perturbation theory in the Appendix, provide 
the framework for studying the issue further. However, 
the relevant (regularized) integrals are unlikely to yield to 
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exact evaluation and thus will have to be studied asymp- 
totically in n, which is beyond the present scope. 



Acknowledgments 

I thank V. Cardoso for correspondence stimulating this 
paper, J. P. Norman for sharing his data, the CUHK 
group for previous collaboration on the Poschl-Teller 
model, and H.Y. Guo for supporting my calculations. 



APPENDIX A: POSCHL-TELLER POTENTIAL 

The PT QNMs are well known 8], but the present 
recapitulation adds a normalization for the eigenfunc- 
tions, as it applies to the perturbation theory which is 
also sketched. Consider 



[dl - V(x) + LJ 2 ]ip(x,uj) = 



V(x) 



V 



df(x/2) 



(Al) 
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and define Q = y 4Vo — j. The outgoing wave to the left, 
f(x—>—oo) ~ 1 • e~ luJX , reads 

= [£(1-0]-*" 

x F{\+iq-2iu, \-iq-2iu; l-2iw; f ) , (A3) 

where F is the hypergeometric function and £ = 
1/ (l+er x ). The outgoing wave to the right, g(x— >oo) 
~ 1 ■ e tuJX , follows as g(x,co) — f{—x,w), i.e., £ i— > 1 — £ 
in (| A3|) . Using the connection formulas for F yields 



g(x,uj) 



T{l-2iu))T{-2iuj) 
T(\+iq-2iu)T(\-iq-2iuj) 

sin[7r(|+ig)] 

-f(x,u>) . 



f{x,-u) 



sin[27riw] 

The Wronskian J(u>) = gf — fg' becomes 

T(l-2iu) 2 



J(w) 



T{\+iq-2iuj)T{\-iq-2iuj) 



(A4) 
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and has zeros 2uj n = ±q — (n+h)i, n = 0, 1, 2, . . . Here, 
we focus on real q and, without loss of generality, on the 
QNMs in the fourth quadrant (plus sign), as has tacitly 
been done in the main text. The derivative is 

J'K) = (-)"+*2m! ^ { _ n _ 2 g (A6) 

= -2w n «/ n5n » , (A7) 

where the second line is a general relation between the 
Wronskian and the norm of the eigenfunctions F° r 
potentials decaying as a series of exponentials [like (|A2|) ]. 
this norm <Ca^> is defined as the constant term in the 



asymptotic expansion of J^dxa(x) for — Xi,X2 — > oo. 
Since the parity of V presently implies that g n = (— )"/„, 



= 2n! 



r(^—n—iq)T(—^—n—iq) 



(A8) 



r(-n-2ig) 

/„(x) = [£(!-£)] "^(-n, -n-2iq; £-n-*g;£) . (A9) 



For the lowest few n, one can verify (| A8|) by substitut- 
ing i|A9() : in terms of £, the integral reduces to analyti- 
cally continued Beta functions [l2j . 

Let us turn to the perturbation theory of (|A1(I . in a 
form which assumes neither completeness of the QNMs 
nor finitely supported perturbing potentials. Since this 
only slightly generalizes previous analyses 0, the 
discussion will be brief. Let V i— > V + eV' + e 2 V" and 
uj n i— > w n + eujj + e 2 w". The starting point is the deriva- 
tive of d), [4 - V + u 2 n \f n = [V[ - 2u n w' n ]f n . The 
solution f n is arbitrary up to two independent homo- 
geneous solutions, say /„ and Xn- The former, corre- 
sponding to an cj-dependent change in normalization, is 
irrelevant; the latter has to be adjusted so that f' n is out- 
going. If the potentials decay at least exponentially and 
Rew„ ^ 0, imposing the boundary conditions does not 
involve any subtleties [l3L Il4| . Satisfying the two condi- 
tions (for x — > ±oo) simultaneously is possible iff 



2w n </2» 



(A10) 



note that the whole calculation is manifestly finite, and 
that neither the 'matrix element' in the numerator nor 
the norm involves complex conjugation [Tll |. Analo- 
gously, one obtains in second order 



<^.f n [{v'-2Lo n uj' n )r n + v"u\> w 



2w„</2> 



2w„ 



(All) 



which may be simpler than its counterpart in 
terms of logarithmic-derivative wave functions |12| . 
One can of course write <^f n (V'—2u! n ui' n )f^^> = 
<^f' n [d 2 —V+Lo 2 ,]f' n ^>, but this does not seem to have an 
advantage. For V even and V odd, the situation of the 
main text, u>' n = 0, but the outgoing condition on f' n 
remains non-trivial. 

A few toy examples, in fact relevant to |JBJ, are 



Vi(x,e) 



V 2 (x,e) 



V 3 (x,e) 



Vn 



ch z (x/2) 

= Vq 
ch 2 [x/(2+e)} 

. 1 2n+l 

>2,n - 



v{ = chr 
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(A12) 



16q 

V 



8 



i : 



UJ 



3 . n 



cK[{x-t)/2] 
. 



V' = 



vl = 



V xsh(x/2) 
~2 ch 3 (x/2) 



U sh(a;/2) 
ch 3 (x/2) 



(A13) 



(A14) 
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For the first few n, w^ 2 ) n can a S a i n be found explicitly 
from IjAlOjl . while uj' 3n — is obvious since f%Vl is odd. 
Working out these examples to the second order <|A11|) . 



and/or studying them asymptotically in n, should be a 
useful (though not yet sufficient) stepping stone for the 
analogous analysis of (|S|)- (fTTJ) . 
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